SPECTRAL THEORY FOR SYMMETRIC ONE-DIMENSIONAL LEVY 
PROCESSES KILLED UPON HITTING THE ORIGIN 

MATEUSZ KWASNICKI 

Abstract. Spectral theory for transition operators of one-dimensional symmetric Levy 
process killed upon hitting the origin is studied. Under very mild assumptions, an 
integral-type formula for eigenfunctions is obtained, and eigenfunction expansion of 
transition operators and the generator is proved. As an application, and the primary 
motivation, integral fomulae for the transition density and the distribution of the hitting 
time of the origin are given in terms of the eigenfunctions. 



1. Introduction 

In two recent articles [T^[T3] . spectral theory for some symmetric Levy processes killed 
upon leaving the half-line was developed. One of the main motivations for these research 
came from fluctuation theory for Levy processes: the distribution of the supremum func- 
tional and the first passage time can be expressed in terms of the eigenfunctions of the 
corresponding transition semigroup. The purpose of the present paper is to obtain similar 
results for processes killed upon hitting the origin, and apply them to the study of the 
hitting time of a single point. The following theorem is our main result. 

Theorem 1.1. Let Xt he a symmetric one- dimensional Levy process, starting at 0, with 
Levy-Khintchine exponent ^(0? '^i^'d suppose that > and 2^\l/"(,^) < for 

^ > 0, and that 1/(1 -|- ^I^(O) ^■^ integrable. Let be the first hitting time of {x}. Then 

^. . 1 pcos^Ae-**W^'(A)FA(x) ,^ , . 

Pt<r, <oo=-/ d\ 1.1 

for t > and almost all x G R. Here Fx is a bounded, continuous function, defined by 

Fx{x) = sm{\Xx\+^x)-Gx{x) (1.2) 

for X G R, where 

and Gx is an ^^^(R) H Co(R) function with (integrable) Fourier transfi 



orm 



^G.K) = cos..(^^-^^^l^) (1.4) 



/or{eR\{-A,A}. 



Remark 1.2. It is easy to check that formula (1.1) can be differentiated in the time 
variable t under the integral sign. By combining it with estimates of Fx, in many cases 
one can obtain asymptotic estimates of the density function of r^., as well as its derivatives, 
in a similar manner as in [13j for the first passage time. This topic will be addressed in 
a separate article. 
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Remark 1.3. When the Levy measure of Xt has completely monotone density on (0, oo), 
then G\{x) is completely monotone on [0, oo) and in many cases can be given by a more 
straightforward formula; see Theorem |1.8[ 



Let us discuss shortly the assumptions of Theorem |1.1[ Since the process Xt is assumed 
to be symmetric, is a real-valued function and ^1/(0 ^ 0. The assumption "^'{C.) > 
and 2^^"(0 < ^'(0 for ^ > is equivalent to the condition V^'IO > 0, < for 

^ > for the function = \I'(v^). This is clearly satisfied by all subordinate Brownian 
motions (and hence for symmetric stable processes and mixtures of such), but also for 
many less regular processes, such as truncated symmetric stable processes. Examples are 
discussed in Section |5} Integrability of 1/(1 + \E'(^)) asserts that the process Xt hits single 
points with positive probability. 

The class of Levy processes studied in |12l [13] consisted of symmetric processes, whose 
Levy measure admits a completely monotone density function on (0, oo). This regularity 
assumption was needed for an application of the Wiener-Hopf method for solving a certain 
integral equation in half-line. For the case of hitting a single point, considered below, 
a more direct approach is available, and therefore much more general processes can be 
dealt with. 

From now on we consider the Levy process Xt starting at a fixed point x G R, and 
denote the corresponding probability and expectation by and E^.. The functions Fx 
in Theorem |1.1| are eigenfunctions of transition operators of Xt killed upon hitting {0}. 
These operators are defined by the formula 



Pf^'^'fix) = E.(/(X,)1 



{t<To}j 



for t > 0, X G R \ 0, and they act on L^(R\0) for arbitrary p G [l,oo]. By A-r\q 

and T>{A-£i\o', L'^) we denote the generator of the transition semigroup P^^'^ acting on 
L^(R \ 0), and its domain; a more detailed discussion of these notions is given in Prelim- 
inaries. Our main results about F\ are contained in the three theorems stated below. 

Theorem 1.4. Let Xt he a symmetric one- dimensional Levy process with Levy- 
Khintchine exponent \E'(0, and suppose that \E''(0 > for^ > 0, and that 1/(1 + ^^(0) 
is integrable. Fix A > 0, and let Fx be defined as in Theorem \1.1\ Then Fx{0) = 0, 
FxeV{AB.\o;L^), and 

An\oFx = -^{X)Fx, P^^^^Fx = e-'^^'^Fx (1.5) 

for t > 0. In addition, for ^ > 0, 



(1.6) 



r ^ . , , ^esin^A + Acos^ 

J oo " WTx' — 

_2cos?9a ^ 2A _ ^'(A) 

~ Jo ¥+e ~ ^(A) - ^(0 

In the sense of Schwartz distributions, 

J^FxiO = cos^Apv ^^^^^^^^^ + TT sin ^,(5,(0 + S^xiO) (1-7) 

(pv stands for the principal value), and 

AFx{x) = -^{\)Fx{x) + ^'{X)cos^xSo{x), (1.8) 
where is the Dirac delta distribution at ^. 



SYMMETRIC LEVY PROCESSES KILLED UPON HITTING 



3 



Remark 1.5. The assertions of the theorem can be interpreted correctly also when 
1/(1 + ^(^)) is not integrable. In this case 'd\ = arctan(oo) = 7r/2, and therefore cos'i^A = 
0, Gx = 0, and finally Fx{x) = cos(Ax), as for the free process. This reflects the fact that 
Xf does not hit from almost all starting points, and therefore P/^^°/(a;) = Ptf{x) for 
almost all a; G R {Pt is the transition semigroup of the free process Xf). 

Remark 1.6. For a given A > 0, the statement of the theorem remains true if the 
condition \1>'(^) > for all ^ > is replaced by the condition \1>'(A) > and \1>(^) 7^ \1>(A) 
for all ^ 7^ ±A. It is an interesting open problem to study the remaining cases, namely: 
when \1/ is increasing on (0, 00) and \E''(A) = 0, and when \l/(^) = ^'(A) has more than 
one solution on (0, 00). 

For estimates and more detailed properties of the eigenfunctions, further regularity of 
the Levy-Khintchine exponent "^l^) is needed. The required assumption is the same as 
in fr2\ [13], and it can be put in the following three equivalent forms. The notion of a 
complete Bernstein function is discussed in Preliminaries. 

Proposition 1.7 (Proposition 2.14 in [I2])- Let Xt he a one- dimensional Levy process. 
The following conditions are equivalent: 

(a) \l/(^) = ipi^"^) for a complete Bernstein function 

(h ) Xt is a subordinate Brownian motion ( that is, Xt = Bz^ , where Bg is the Brown- 
ian motion, Zt is a non- decreasing Levy process, and Bg and Zf are independent 
processes), and the Levy measure of the underlying subordinator Zt has completely 
monotone density function; 

(c) Xt is symmetric and the Levy measure of Xt has a completely monotone density 
function on (0,oo). 



00 



00 



Theorem 1.8. Let A > 0. Suppose that the assumptions of Theorem I.4 hold true. 

(a) //2^^"(0 < ^'(0 for all ^ > 0, then we have i^x G [0,vr/2) and \Gx{x)\ < 
Gx{0) = sm^x- 

(h) Suppose that any of the equivalent conditions of Proposition 1.1 holds true. Then 
■i^A G [0,vr/2), Gx is completely monotone on [0,oo), andGx{x) = £7a(|x|) (where 
C is the Laplace transform) for a finite measure 7a. In particular, Gx{x) > 0, 
and Gx is integrable, 

G.Wdx = _FG.(0) = ^^(l-A|^). (1.9) 

// in addition \1/ extends to a function ^E'^(0 holomorphic in the right complex 
half-plane {.^ G C : Re^ > 0} and continuous in G C : Re.^ > 0}, and 
^+(i^) ^ $(A) for all ^ > 0, then for all x eR, 

^^^"^^ TT I vI/(A)-M/+(20 ^' ^ ' 

Note that the assumptions for part (a) of the theorem are exactly the same as in 
Theorem II. 1[ 

Before the statement of the third main result about the eigenfunctions Fa, let us make 
the following remark. The 'sine-cosine' Fourier transform J-", defined by the formula 

(poo POO 
/ f{x) cos{Xx)dx, / f{x) s'm{Xx)dx 
J —00 J —00 

for / G Cc(R) and A > 0, extends to a unitary (up to a constant factor ^/^^) mapping 
of i^^(R) onto L^((0,oo)) © L^((0, 00)), in which the free transition operators Pt take a 
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diagonal form. Namely, Pt transforms to a multiplication operator e"***^'^-', that is, 

:FPJ{X) = e-**W J-/(A); 



see Remark 2/7 The above formula is an explicit spectral decomposition, or generalised 
eigenfunction expansion of Pf. The next theorem provides a similar result for transition 
operators P^'^'^ of the killed process. The corresponding transform is given by a similar 
formula as J-", with cos(Aa;) replaced by Fx{x). 



Theorem 1.9. Suppose that the assumptions of Theorem I.4 are satisfied, and let Fx be 
the eigenf unctions from Theorem \1.4\ Define 



neven/(A) = / f{x)Fx{x)dx, 



—00 
00 



nodd/(A) = / f{x) sm{Xx)dx, 
J —00 

n/(A) = (neven/(A),nodd/(A)) 

for f G Cc(R \ 0) and A > 0. Furthermore, define 

POO POO 

n*(/i,/2)(x)= / /i(A)F,(x)rfA+ / /2(A)sin(Ax)rfA 
JO Jo 

for /i, /2 G Cc((0, oo)) and a; > 0. Then tt^^^'^H and tt^^^"^]!* extend to unitary operators, 
which map L^{R\0) onto L^{{0,oc)) ® L'^{{0,oo)), and L^{{0,oo)) ® L'^{{0,oo)) onto 
L^(R\ 0), respectively. Furthermore, UIl* = tt, 11*11 = tt, and 

npf \V(A) = e-*'^wn/(A), (1.11) 

for all / G L^(R \ 0). In addition, f G L^(R \ 0) belongs to X'(^r\o; L'^) if and only if 
\l/(A)n/(A) is square integrable on (0, oo), and in this case for A > 0, 

n^R\o/(A) = -vI/(A)n/(A). (1.12) 

In a rather standard way, explicit spectral representation for P^ yields a formula for 
the kernel of the operator, that is, for the transition density. 

Corollary 1.10. Suppose that the assumptions of Theorem \1.4\ are satisfied, and 
2^\l/"(,^) < "^'l^) for all C, > 0. Then the transition density of the process killed upon 
hitting the origin (that is, the kernel function of P^"^^) is given by 

pf\°(x,i/) = - / e-'''^'\Fxix)Fxiy) + sin(Ax) sin(Ai/))rfA 

Jo 

+ _ / e-'''^'^Fxix)Fxiy)dX. 

Jo 



pt{x - y) - pt{x + y) 1 i^-^^) 



2 

Formula (1.13), however, is of rather limited applications due to many cancellations 
in the integral with two oscillatory terms. Nevertheless, it is one of the very few explicit 
descriptions of the transition density of a killed Levy process. 

By symmetry, we have PxiTo > t) = Poi^x > t). Hence, Theorem |1.1| follows trivially 
from the following result. In the theorem, an extra condition on ^' is required in order 
to use Theorem ll.Slfa). 



Theorem 1.11. Suppose that the assumptions of Theorem L4_ are satisfied, and 
2C*"(0 < ^'(0 for all^>0. Then 

,1 /""^ cos^9Ae-**(^)^'(A)FA(x) „ 
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for t > and almost a// x G R \ 0. 

Since Pxi'^o > t) = Jii\QpT^'^{^^y)d-y^ oii^ could naively try to derive formula (1.14) 



by integrating (1.13) over y G R \ 0. Heuristically, changing the order of integration 



and substituting J-'Fa(O) = cos 'i9A^'(A)/^(A) for the integral F\{x)dx (which is not 



absolutely convergent, or even convergent in the usual way) would yield (1.14) with 
Px{t < To) instead of Px{t < Tq < oo). Hence, during this illegitimate use of Fubini, 
Po{to = C)o) is lost. 

Relatively little is known about the distribution of the hitting time of single points, r^. 
In the general case, the double integral transform (Fourier in space, Laplace in time) is 
well-known. 



/r ^ + ^(0 

for 2; > and ^ G R, where c{z) is the normalisation constant, chosen so that the integral 
of the right-hand side is 1 (see Preliminaries). Hence, 

uJx) 



ZTO 



for z > and x G R, where Uz{x) is the resolvent (or ^-potential) kernel for Xt. However, 
Uz{x) is typically given by an oscillatory integral (namely, the inverse Fourier transform of 
l/(2-|-\I'(^)), and therefore it is not easy to invert the Laplace transform in time, even for 
symmetric a-stable processes. Some asymptotic analysis of Px{to > t) (and many other 
related objects) can be found, for example, in |Hl [HI ED]. Some more detailed results in 
this area have been obtained for spectrally negative Levy processes, see, for example, [6] 
and Section 46 in |16) . For totally asymmetric a-stable processes, a series representation 
for Px(to > t) was obtained in jl5] . 

The hitting time of a single point is closely related to questions about local time and 
excursions of a Levy process away from the origin, see [5l [T9j. Similar questions arise 
when killing a Levy process is replaced by penalizing it whenever it hits 0, cf. |2T]. In 
fact, using the same method, one can find generalized eigenfunctions for the transition 
operators of the penalized semigroup. Finally, spectral theory for the half-line has been 
succesfuUy applied in the study of processes killed upon leaving the interval |10| [TTl [T3] 
(see also [ll!2j) and higher-dimensional domains [8j. One can expect similar applications 
of the spectral theory for R \ 0, developed in the present article. 

The remaining part of the article is divided into four sections. In Preliminaries, we 
describe the background on Schwartz distributions, complete Bernstein and Stieltjes func- 
tions, and Levy processes and their generators. Section |3] contains the derivation of the 
formula for the eigenfunctions Fx- The generalised eigenfunction expansion and formulae 
for the transition density and the distribution of the hitting time are proved in Section |4j 
Finally, examples are studied in Section [5] 

2. PRELIMINARIES 

2.1. Schwartz distributions. Some background on the theory of Schwartz distributions 



is required for the proof of Theorem 1.4 A general reference here is [H]; a closely related 
setting is described with more details in |12] . 

The class of Schwartz functions in R is denoted by S, and S' is the space od tempered 
distributions in R. If G 5 and F G S', we write {F,(p) for the pairing of F and ip. 
The Fourier transform of (/9 G 5 is defined by J^f{C,) = e'^^^ip{x)dx. For F G S', J^F 
is the tempered distribution satisfying {J^F, ip) = {F, J^ip) for all ip & S. This definition 
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extends the usual definition of the Fourier transform of ^^^(R) or ^^^(R) functions and 
finite signed measures on R. 

The convolution of (fi, ip2 E S is defined as usual by ipi*ip2{x) = ipi{y)ip2{x — y)dy. 
When if & S and F e S', then F * (f is the infinitely smooth function defined by 

F * Lp{x) = {F,Lp^), where ip^{y) = ip{x - y). 

The convolution does not extend continuously to entire space S'. We say that -Fi, F2 G S' 
are S'-convolvable if for all cpi, ip2 G S, the functions Fi * (pi and F2 * ip2 are convolvable 
in the usual sense, i.e. the integral J_^{Fi * ipi){y){F2 * f2){x — y)dy exists for all x. 
When this is the case, the S' -convolution Fi * F2 is the unique distribution F satisfying 
F*ipi*ip2 = {Fi* (pi) * {F2* ip2) for 'Pi,(p2 £ <S- There are other non-equivalent definitions 
of the convolution of distributions; see [18j. 

The support of a distribution F is the smallest closed set supp F such that {F, ip) = 
for all (/? G iS such that ip{x) = for x G suppF. If any of the tempered distributions 
Fi, F2 has compact support, or if Fi * (f is bounded and F2 * ip2 is integrable for all 
<y9 G iS, then Fi and F2 are automatically iS'-convolvable. Note that ^'-convolution is not 
associative in general. 

If distributions Fi, F2 are iS'-convolvable, then ^-'(-Fi * F2) = J^Fi ■ TF2 (the exchange 
formula), where the multiplication of distributions extends standard multiplication of 
functions in an appropriate manner. Below only a few special cases are discussed, and 
for the general notion of multiplication of distributions, we refer the reader to |18| . 

When Fi is a measure and F2 is a function, then their multiplication Fi ■ F2 is the 
measure F2{x)Fi{dx), as expected. Next example requires the following definition. By 
writing Fi = pv(l/(a; — xq)) we mean that 

(Fi,(^)=pv/ ' dx = lim / dx; 

the principal value integral pvj is defined by the right-hand side. When Fi = pv(l/(a; — 
Xq)) and F2 is a function vanishing at Xq and differentiable at Xq, then Fi-F2 is the ordinary 
function F2{x) / {x — Xq) . The third example is the extension of the previous one. Let / be 
a function vanishing at a finite number of points Xi, X2, x„, continuously differentiable 
in a neighborhood of each Xj with f'{xj) 7^ 0, and such that is bounded outside 

any neighborhood of {xi, In this case, Fi = pv(l//(x)) is formally defined by 

(Fi, ^) =pvr^dx = hm / ^ dx. 

Clearly, Fi = Yl^=i f'i^j) Pv(l/(^ ~ Xj)) + F, where F is a genuine function. Hence, by 
the previous example, we obtain the following result: when Fi = pv(l//(x)) is as above 
and F2 is a function which vanishes at each Xj and is differentiable at each Xj, then Fi ■ F2 
is the function F2{x)/f{x). 

The Laplace transform of a function F on [0, 00) is defined by CF{^) = F{x)e~^^dx. 
For a (possibly signed) Radon measure m on [0, 00), we have Cm[S,) = e~^^m{dx). 



2.2. Stieltjes functions and complete Bernstein functions. A general reference 
here is [T7]. A function f{x) is said to be a complete Bernstein function (CBF in short) 
if 

f{x) = c^ + C2X + - — x>0, (2.1) 

TT /n X + s s 
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where c\ > 0, C2 > 0, and m is a Radon measure on (0, 00) such that 
J min(s~^, s~^)m((is) < 00. The right-hand side clearly extends to a holomorphic func- 
tion of X G C \ (— 00,0], and we often identify / with its holomorphic extension. A 
function /(x) is a Stieltjes function if 

Ci 1 /""^ 1 
f{x) = - + C2 + - m{ds), x>0, (2.2) 

X IT Jq X + S 

where Ci > 0, C2 > 0, and m is a Radon measure on (0, 00) such that J min(l, s~^)m{ds) < 
00. Below we collect the properties of Stieltjes and complete Bernstein functions used in 
the article. 

Proposition 2.1 (see |17)). (a) A function f is a CBF if and only if f{z) > for z > 
and either f is constant, or f extends to a holomorphic function in C \ (—00, 0], 
which leaves the upper and the lower complex half-planes invariant. 

(b) A function f is a Stieltjes function if and only if f{z) > for z > and either f 
is constant, or f extends to a holomorphic function in C\ (— oo,0], which swaps 
the upper and the lower complex half-planes. 

(c) If f is positive, then the following four conditions are mutually equivalent: f{x) 
is a CBF, x/f{x) is a CBF, l/f{x) is a Stieltjes function, f{x)/x is a Stieltjes 
function. 

(d) A function f is a Stielties function with representation ( 2.2[ ) if and only if f is 
the Laplace transform of the measure C2SQ{ds) + (ci + Cm{s))ds. □ 

Proposition 2.2 (Proposition 7.14(a) in [17j). If ip is a non-constant CBF, then 

MO - , (2.3) 

(extended continuously at A^, so that ^/'a(A^) = 4j{X'^)/{X'^ip'{X'^))) is again a CBF. □ 

Proposition 2.3 (Corollary 6.3 in ^7J, Proposition 2.18 in [12j). 
(a) Let f be a CBF with representation (2.1). Then 

ci = lim f{z), 02 = lim ^ , (2.4) 

z-^0+ 2->oo Z 



and 



mids) = lim (Im /(— s + z£:)(is), (2.5) 



with the limit understood in the sense of weak convergence of measures, 
(b) Let f be a Stieltjes function with representation (2.2). Then 



ci = lim {zf{z)), C2 = lim f{z), (2.6) 

2-S.0+ 2^-00 



and 



m{ds) + nc2So{ds) = lim {—Im f{—s -\- ie)ds), (2.7) 

e-s>0+ 

with the limit understood in the sense of weak convergence of measures. □ 

Proposition 2.4 (Proposition 2.21 in [I2j)- (d) If f is a complete Bernstein func- 
tion, then \f{z)\ < C{e)f{\z\) and \f{z)\ < C{f,e){l + \z\) for z E C, 
I Aigzl < IT — e, e e (0,7r); 
(b) Iff IS a Stieltjes function, then \f{z)\ < C{e)f{\z\) and \f{z)\ < C{f,e){l + \z\-^) 
for z e C, \ Argz\ < TT — e, e E {0,Tc). □ 
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2.3. Levy processes. Below we recall some standard notions related to Levy processes 
and describe the setting for the present article. General references for Levy processes 
are P, |16|; for the properties of semigroups of killed Levy processes, see fl2j and the 
references therein. 

Throughout this article, Xf is a one-dimensional Levy process with Levy-Khintchine 
exponent ^(0; that is, 

Eoe'«^* = e-**(«), t > 0, e e R. (2.8) 

The probability an expectation corresponding to the process starting at a; G R are denoted 
by Px and Ej,, respectively. The following assumptions are in force throughout the article: 

(A) Xt is a symmetric process; 

(B) Xf hits every single point in finite time with positive probability. 

The first assumption is equivalent to the condition "^l^) = ^(—0 ^ ^"^^ 

to the condition and > for ^ G R. For symmetric processes. Assumption (|B]) is 

equivalent to integrability of l/{z + ^I'(O) over ^ G R for some (and hence for all) ^ > 0, 
see Theorem 2 in [9j, Theorem 11.19 in [3j or Theorem 43.3 in fT6]. This equivalent form 
is more convenient in applications. 

Both conditions are natural for the spectral theory. When Assumption |A] fails, i.e. 
when Xt is not symmetric, the operators P^"^^" are not self-adjoint (in fact, they even fail 
to be normal operators). Assumption B asserts that the operators P^"^^ are not equal to 



Pf] the spectral theory for the latter is trivial, see Remark 2.7 



The transition operators Pt of the free process X^ are defined by 

PJ{x) = EJiXt) = [ /(y)P.(Xi G dy), t>0,xeR, 



JR 

whenever the integral is absolutely convergent. The operators Pt (acting on L^(R)) are 
Fourier multipliers with Fourier symbol e"**'-^-*, and hence they are convolution operators. 
By assumption ([B]) and Fourier inversion formula, the convolution kernel is an integrable 
function pt{x), called the transition density. 

Let D C R be an open set, and let r^) = infjt > : Xt ^ be the first time the 
process Xt exits D. The transition operators of the process Xt killed upon leaving D are 
given by 

Pt^'fix) = E,{f{Xt)lt<r,) = [ f{y)P.iXt edy;t<Tn), t>0,xeD. 

The kernel function pj^{x,y) of (the transition density of the killed process) exists, 
and < (x, y) < pt{x — y) {t > 0, x,y ^ D). 

By L^{D), p G [1, oo), we denote the space of real-valued functions / with finite p-norm 
||/||lp(D) = (/f) \ f{x)\^dxy^^. The space of essentially bounded functions on D is denoted 
by L°°{D). The space of bounded continuous functions is denoted by Cf,(R), and Cq(R) 
is the set of continuous functions vanishing at ±oo. By Cq{D) we denote the space of 
Co(R) functions vanishing on the complement of D. 

In this article, we only consider D = R \ 0, so that To = tq is the hitting time of the 
origin. In this case clearly Lp(R\ 0) can be identified with L''^(R); nevertheless, we use 
the notation Lp(R \ 0) to emphasise that the corresponding operators are related to the 
process killed upon hitting the origin. 

The following properties of transition semigroups are described in detail in [12]. The 
operators Pj and P^"^^" form a contraction semigroup on L^{IV) and Lp(R\0), respec- 
tively, where p G [1, oo]. This semigroup is strongly continuous when p < oo. Also, Pt is 
a contraction semigroup on C;,(R) and Co(R), strongly continuous on the latter space. 



SYMMETRIC LEVY PROCESSES KILLED UPON HITTING 



9 



By Theorem 11.19 in [3j or Theorem 43.3 in |16] . is regular for itself, and so p^^^ is a 
strongly continuous semigroup of operators on Co(R \ 0) (see also |12]). 

We denote the generators of the semigroups Pt and P^ (on appropriate function 
spaces) by A and ^r\o, respectively. The corresponding domains are denoted by T>{A; X) 
and 'D{A^\q\ X), respectively, where X indicates the underlying function space. For ex- 

R,\ 

ample, 'D{Ayi\q'-, L°°) is the domain of the generator of the semigroup of contractions 
on L°°(R \ 0), that is, the set of those / G L°°(R \ 0) for which the limit 

An\of = lim ' ; ^ (2.9) 



exists in L°°(R \ 0). Note that the limit (2.9), if it exists, does not depend on the choice 



of the underlying function space. Therefore, using a single symbol A for operators acting 
on different domains T>{A; X) causes no confusion (again, see [121 ) . 

The resolvent kernel, or the ^-potential kernel, is defined by Uz{x) = pt{x)e^^^dt. 
For each z > 0, it is a bounded, continuous function of a; G R, and J^Uz{^) = l/(-2 + ^(.^)). 
Furthermore, Uz{x) is the convolution kernel of the inverse operator of ^ — .4 (acting on 
L^(R)). On the other hand, by Theorem 11.19 in [3j or Theorem 43.3 in fTU\, Uz{x)/uz{0) 
is the Laplace transform of the distribution of tq, 

E.e— « = ^, z>0,xER. 

Uz{0) 



This identity will be used in the proof of Theorem 1.11 

Following [12J, we introduce the distributional generator of Xt, denoted by an italic 
letter A: we let A be the tempered distribution satisfying Aif = A * (p for all (f E S. 
By symmetry of Xf, equivalently we have {A,(f) = Aip{0) for (f E S. The distributional 
generator can be thought of as a pointwise extension of the generator A. The following 
result from |12| describes the connection between A and ^r\o, the generator of the killed 
semigroup. 

Lemma 2.5 (see Lemma 2.10 in [12]). Let F G Co(R\0). // the distribution A* F is 
equal to a Co(R\ 0) function plus a distribution supported at {0}, then F G T>{Aii\o; C*o) 
and An\oFix) = A* F{x) for X eR\0. □ 

The following technical result has been proved in |12| for general domains. We record 
its version for R \ 0. 

Lemma 2.6 (see Lemma 2.11 in [12]). Let F G C^(R), that is, F an all derivatives of 
F are bounded functions on R. Suppose that for some r > 0, F{x) = for x G [— r, r], 
and AF{0) = 0. Then F G V{Ar\o; Co) and An.\oF{x) = AF{x) for all x eR\0. □ 

Remark 2.7. The spectral theory for the free process is very simple, thanks to the Levy 



Khintchine formula (2.8). Indeed, the function e*'^^ (A G R) is the eigenfunction of Pt and 
A, with eigenvalue e"*^'-'*''' and — ^(A), respectively. Clearly, these eigenfunctions yield 
the generalised eigenfunction expansion of Pf and A by means of Fourier transform. 

Since Xt is symmetric, ^I'(A) = \l/(— A) is real. It follows that sin(Ax) and cos(Ax) 
(A > 0) are also eigenfunctions of Pt and A. Consequently, sin(Ax) yields generalised 
eigenfunction expansion for odd L'^(R.) functions (via Fourier sine transform), while 
cos(Ax) plays a similar role for even -^^^(R) functions (via Fourier cosine transform). □ 

3. Eigenfunctions in R \ 
In this section the formula for Fx is derived, and some properties of eigenfunctions are 



studied. More precisely, we prove Theorems 1.4 and 1.8 The argument follows closely 
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the approach of [12], where the case of the half-hne (0, oo) is studied. Noteworthy, in 
our case there is no need to employ the Wiener-Hopf method, which makes the proof 
significantly shorter and simpler. 

Mimicking the definition of distributional eigenfunctions in half-line (Definition 4.1 
in [12]), we introduce the notion of distributional eigenfunctions of ^r\o- Note that the 
condition F[0) = has no meaning for general Schwartz distributions F, so in contrast to 
the definition of |12] . at this stage we do not assume that F vanishes in the complement 
of the domain (that is, at the origin). 

Definition 3.1. A tempered distribution F is said to be a distributional eigenfunction 
of ^R\o, corresponding to the eigenvalue — \E'(A), if F is 5'-convolvable with A, and 
A* F + 'ip{\^)F is supported in {0}. 

Note that, in particular, cos(Aa;) is a distributional eigenfunction of v4r\o- However, 
it is not the one we are looking for, as it does not vanish at 0. By copying the proof of 
Lemma 4.2 in [T2| nearly verbatim, one obtains the following result. We only sketch the 
prove and omit the technical details, referring the interested reader to [12j . 

Lemma 3.2. Let A > 0, and suppose that F is a distributional eigenfunction of A-^\q, 
corresponding to the eigenvalue — \E'(A). If F is bounded, continuous on R, -F(O) = and 

lim (F(x) -Csin(|Aa;| +^9)) = 

s— >-±oo 

for some C,d eH, then F G V{A^\q\ L°°) and A^\qF = -^(A^)^. 

Sketch of the proof. It is possible to find an infinitely smooth function /i such that 
/i(x) = when |x| < 1, /i(x) = Csin(|Ax| + when |x| is large enough, and 



Afi{0) = 0. We let /2(x) = F{x) - fi{x). By Lemma [2^ /i e V{An\o;Cb) and 
'^R\o/i(2;) = -^fii^) = ^ * fi{x)- With a little effort, one shows that Lemma 2.5 
applies to /2. It follows that /2 G ^^(^r\o; C'o) and Aii\of2{x) = A* f2{x). Hence, 
F = f^ + f^e V{Ar\o; Cb) and Ar\oF = A* F = -ip{\^)F, as desired □ 

It is relatively easy to find a formula for distributional eigenfunctions of ^r\o satisfying 



the assumptions of Lemma [3. 2 [ First, we give a brief, heuristic derivation of the formula. 
Suppose that Fx is a bounded, continuous, even function on R such that F\{0) = and 
A * F + %1){\^)F is supported in {0}. A tempered distribution is supported in {0} if 
and only if its Fourier transform is a polynomial. Hence, (— \E'(^) + '^{X))J^F\{E,) is a 
polynomial Q{C)- It follows that the distribution J^Fx is expected to have form 

^{X) - vi/(o 

(the principal value corresponds to singularities at ±A), plus some distribution supported 
in {—A, A} (the zeros of \E'(A) — \E'(^)). The function Fx should be as regular as possible, so 
we assume that Q is constant (say, Q{^) = aA^I/'(A)), and that the distribution supported 
in {—A, A} is a combination of Dirac measures (say, 7ibx{Sx + This suggests the 

following definition: 

TFxiO = axpy ^ ""^'^^'^^^ + ^-'^^^^ 

for some ax,bx G R. We can normalize this definition by assuming that a\ + b\ = 1 and 
o-x > 0, so that ax = cos'dx and 6a = sint^A for some 'dx £ [~^/2,7r/2). Furthermore, 
the condition Fx{0) = can be formally rewritten as J-'Fx{^)d^ = 0, which gives 
a linear equation in ax and 6a- After soving this equation, we obtain formulae given in 
Theorem II. 4[ 
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Remark 3.3. Before the proof Theorem |1.4[ let us make tho following observation. The 
primitive function of 2\/{\^ — ^'^) is log((A + ^)/(A — ^)). This implies that for A > 0, 



pv 



TT 



2A 







d^ = 0. 



By a similar direct calculation (we omit the details), for A,.^ > 0, 



pv 



TT 



2A 







e + A2 - C' 



dc 



A 



A2 + e 



Hence, formulae of Theorem 1.4 take a simpler form 



arctan 



1 



TT 



pv 



and 



FA(x)e-«l^ldx = 2 



^(A)-^(O 
^ sin {}x 



(3.1) 



2 cos dx 
H pv 



e + A2 

oo ^ 



*'(A) 



(3.2) 



e + e^(A)-^(c) 



Proof of Theorem l.^. Since \I^(^) is smooth near ,^ = A, the integrand in (1.3) is a 
bounded function, and the integral is finite by the assumption. By a similar argument, 
the right-hand side of (1.4) is a bounded integrable function. Hence, (1.4) indeed defines 
an i^2(R) n Co(R) function Gx{x), and so Fx{x) is well defined and belongs to Cfe(R). 
Denote 



2 roo 
pVj 



TT 



^>'{\)di 







v^(A) - v[/(0 
so that dx = arctan (/Ca) (by ( |3.1 )). Let 

1 



ax = COS-i?; 



sin-i?; 



A > 0, 



(3.3) 



(3.4) 



Note that J^Gxii) = axi2\/{y - - ^'(A)/(^(A) - ^(0))- We define 

F^(x) = sin(|Aa;| + i!}x) = ax sin |Ax| + bx cos(Ax), 

so that Fxix) = F^{x) - Gxix). 
By Fourier inversion formula. 



GxiO) 



1 

2^ 

OA 
TT 



TGxiOd^ 
2A 



pv 



A2-^' 



d^ - pv 



*'(A) 



^(A)-^(O 



d^ 



The first principal value integral vanishes (cf. Remark 3.3). It follows that 
Gx{0) = axKx = cos{^x) tan(^A) = sin(^9A) = Fa*(0), 

and therefore -Fa(O) = 0, as desired. 

Note that 2tt5x{C) is the Fourier transform of e~*'^^, while pv2i/(A — ^) is the Fourier 
transform of — e~*^^ signx. Hence, tt{5-x{C) + '^a(O) ^^e Fourier transform of cos(Ax), 
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and pv(2A/(A^ — ^^)) is the Fourier transform of sin |Ax|. It follows that, 



7r6A(5-A(0+^A(0)+«APV 



2A 



OA 



2A 



ax pv 



A2-e ^(A)-vI/(0 
+ vr6A(5A(0 + 5-A(0), 



and (1.7) is proved. As any bounded function, Fa is iS'-convolvable with A, and by the 
exchange formula we have (see Preliminaries) 

J^{A * Fa + vI/(A)Fa)(0 = (-^(0 + ^(A))^Fa(0 = ax^'{\). 



In particular, A * Fx + \E'(A)Fa = aA^''(A)5o; which proves (1.8). Formula (3.2) follows 
from (1.7) by the exchange formula, since the Fourier transform of e""^'^' is 2a/ [a? + ^^). 
By Remark 3.3, (1.6) follows. It remains to prove that Fa G V{A^\q] L°°) and that ( 1.5[ ) 
holds true. 

By ( |1.8 ), Fa is a distributional eigenfunction of ^r\o- Furthermore, Fa G Cfc(R 



Fa(0) = and Fa(x) — F^{x) = Gx{x) converges to as x — )■ ±oo. By Lemma 3.2 



Fa G V{An\o;L°°) and .4r\oFa = -^(A)Fa. Finally, for x G R \ and t > we have 
(see, for example, [7]) 



P^\'Fx{x) 



Fx{x) + / P!'\'An\oFx{x)ds 
Jo 

Fx{x)-^i\) [ P]'\^Fxix)ds. 
Jo 



For a fixed x G R \ 0, this integral equation is easily solved, and we obtain P^'^'^Fx{x) = 
e~**^'^^FA(x), as desired. □ 

Suppose that ^1/(0 = '^(^^) fo^' some function (^ > 0). Note that with the notation 
of (2.3 ), we have 



J'GxiO = ax 



2A 



2AV^'(A2 



2A 



ax 



1 - 



^a(A2 



(3.5) 



A (^A)A(e) 



Proof of Theorem 1.8 . Part (a). Let ipi^) = \E'(v^). Then we have > and 

^"(0 — fc)^ > 0; hence, is increasing and concave. First, we prove that 'dx > 0. 
Let Kx be defined by (3.3). Observe that since ipxi^"^) is infinitely smooth near ^ = A, 
we have 



1 AMA) r 

vr vl/(A) 



JO 

1 , A2^'(A) 
lim ^ , 

TT q^l- ^(A) 



qX 



+ 



X/q 



Me) 



d^. 
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By substituting ^ = As for ^ < A and ^ 
integrals, we obtain that 



A/s for ^ > A, and then combining the two 



1 AM/'(A) f^{^p,{xys')-M>^'s')) 



IT 



^(A) 



ds. 



Since ijj is concave, ipx is a non-decreasing function (e.g. by a geometric argument: 
\~'^'4){\^)/iIj\{C) = (V'(C^) ~ ''/'(A^))/(^^ — A^) is the difference quotient of tp, and hence a 
non-increasing function of ^). It follows that Kx > 0, and therefore -^x = arctan(i^'A) ^ 0. 
Monotonicity of tpx and formula ( 3.5[ ) imply also that J-'Gx{^) > for all ^, which, by 
Fourier inversion formula, proves that |G'a(OI — Gx{0). 

Part (b). Suppose now that defined as above is a complete Bernstein function. 
Clearly, ip is concave, so, by part (a), "t^x ^ [0,vr/2). The proof of (1.10) is very similar 
to the proof of formula (1.4) in [12j. 



Let Gx{x) 



L(0,oo) 



'x)Gx{x). We will show that CGx is a Stieltjes function. For a > 0, 



the Fourier transform of e "^'^1 is 2a/ (a^ + ^'^) 
we have 

'C-GxiO = ^ / e'^^^^Gx{x)dx 



Hence, by Plancherel's theorem, for > 
^J'G xiO 



1 



n 



Substituting ( = s^/^, we obtain 



By (3.5), 



CGxiO = 

TT 



J'GxiO 



V^TGxisV^) 



2 cost!}) 



1 



ds. 



dC 



(3.6) 



By Proposition 2.2 
sition 



2.1 



A Mxie) ■ 

ipx and {'ipx)x are complete Bernstein functions. Hence, by Propo- 



1/(V'a)a is a Stieltjes function, bounded on (0,oo), and so J-'Gx extends to 
a holomorphic function in the right complex half-plane Re^ > 0. Abusing the notation, 
we denote this extension with the same symbol J-'Gx- By Proposition |2.4[ b), J-'Gx{^) is 
bounded on every region —it/ 2 + e < Aig^ < 7i /2 — e. It follows that the right-hand 
side of (3.6) defines a holomorphic function in C \ (— oo,0], which is the holomorphic 
extension of CGx- We denote this extension by g{C,)- 

Assume that Im,^ > and Re^ < 0. By substituting s - 



'CV^ (C is a complex 



variable here) in the right-hand side of (3.6), we obtain 

J'Gx 



9(0 



ec) 



dC- 



1 + C 

The integrand has a simple pole at C = i, with the corresponding residue 
{—i/2)J^Gx{—iC)- By an appropriate contour integration, the residue theorem and pas- 
sage to a limit (for the details, see Lemma 3.8 in [12]), 



TC 



J'Gxi-^O 



1 + e 

When Im^ > 0, Re^ > 0, we have in a similar manner 



dC + J'GxHO- 



9(0 



TT 



j'Gxm 



dC- 



(3.7) 



(3.8) 



Recall that -FGa(C) = (SA'^ cos79a)/(V^a)a(C^ 



, and l/{ipx)x is a Stieltjes function. Hence, 



Im J'Ga(C) > when ReC > and ImC > 0, and Im J'G'a(C) < when ReC > and 
Im^ < 0. By Proposition 2.1 b), formulae (3.7) and (3.8) imply that g = CGx is indeed a 
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Stieltjes function. Hence, by Proposition 2.1 d), G\ is the Laplace transform of a Radon 



measure 7^ on [0, 00). Since Gx G L'^(R.) H L°°(R), 7^ does not charge {0}, and it is a 
finite measure. 

The assumption for the last part of theorem can be rephrased as follows: the holomor- 
phic extension of to the upper complex half-plane has continuous boundary limit on 
(—00,0), which will be denoted by and 7^ ^(A^) for all ^ > 0. In this case, 

by ( 3.7[ ) and (3.5), the continuous boundary limit of g on (— oo,0) approached from 



the upper half-plane exists, and it satisfies 

+ 1 r J'GM) , / 2A 2A^'(A2^ 
9 {-0 = / 1 , >o dC + ax' 



By Proposition |2.3Kb), it follows that 



^xidO = - Img i-Od^ = Im ^^^,^_^^^_^,^ dC. 



Formula ( 1.10[ ) is proved. □ 

4. ElGENFUNCTION EXPANSION IN R \ 



In this section we prove Theorem 1.9, which states that the system of generalised 



eigenfunctions Fx{x) and sin(Ax), A > 0, is complete in L^(R\0). Throughout this 



section assumptions of Theorem 1.4 are in force, that is, '^{^) is the Levy-Khintchine 
exponent of a symmetric Levy process Xt, 1/(1 + ^I^(O) is integrable, and ^^'(0 > for 
^ > 0. Our argument is modelled after the proof of Theorem 1.9 in jl3], providing a 
similar result for the half-line. Noteworthy, in contrast to the previous section, here the 
case of R \ appears to require essentially more work than the half-line (0, 00). 

We begin with some auxiliary definitions and four technical lemmas. For z E C \ 
(—00, 0], we define 

1 f°° 1 

^i^,z) = - dC, 
^{^1, 6, z) = - -^777^—. 



The integrals are convergent by integrability of 1/(1 + \l/(^)). By a substitution \I'(C) = s, 
one easily sees that '^{z), ip{C,, z) and v^(^i, ^2, z) are Stieltjes functions of z. Representing 



measures (as in Proposition 2.3 ) for these and related Stieltjes functions play an important 
role in the sequel. 

We remark that the above three functions are related to the resolvent (or z- 
potential) kernel Uz{x) of the operator A., for example, (p{z) = Uz{0) and ^p{C,,z)/C, = 
Uz{x)e~^^^^dx (see Preliminaries). These connections are only used in the proof of 
Theorem 11.111 

Lemma 4.1. For any ^ > 0, the function ip{^^z)/ (p{z) is a Stieltjes function of z. 



Proof. Clearly, Lp{$^, z) /ip{z) > for 2; > 0. By Proposition 2.1[ b), it suffices to show 



that lm{(p{^, z)/(p{z)) < when Imz > 0; then automatically lm{ip{^, z)/(p{z)) 
— lm.{(f{^, z)/(f{z)) > when Imz < 0. Let z = x + iy for real x, y. We have 
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Hence, by expanding the integrals and a short calculation, we obtain 

ip{^, z) Im ip{i, z) Re ip{z) - Re ip{i, z) Im (^{z) 
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Im 



ip{z) 



■K'^\i^{Z) 



\^{Z)\^ 

e 



■X 



e + cf 



((v|/(Cl)+x)2+y2)((vl/((2)+x)2+y2; 



Another version of the above formula is obtained by changing the roles of Ci and ^2- By 
adding the sides of the two versions of the formula, we obtain a symmetrised version. 



2Im 



ip{z) 



e 



e 



^'I'pizWJo Jo \e+cf e+ci 

l/(vI/(Cl)-Vl/(C2)) 



X 



((vl/(Cl)+x)2+y2)((vI/(C,)+a;)2+^2 



dCidC2- 



Note that ^(C) is increasing on (0, oo), while C,'^/{^'^ + C^) decreases with C > 0. It 
follows that when y > 0, the integrand on the right-hand side is non-positive, and so 

lm{^{^,z)/^{z))<0. □ 



Lemma 4.2. For any ^i,C,2 > 0, the function 



viz) 



is a Stieltjes function of z. 



Proof. The argument is very similar to the proof of Lemma For simplicity, we denote 
z) = ai + ibi, ip{^2, z) = a2 + ib2, f{z) = a + ib for real a, ai, 02, b, bi, 62- We have 



Im 



'fi^i, z)^p{^2, z) (0162 + 6102)0 - (aia2 - bib2)b 



ip{z) 



a2 + 62 



(a^ + 6^)6162 — {obi — aib){ab2 — a2b) 
6(a2 + 62) 

lYaip{ii,z) Imv?(^2,^) ^ Im(v9(6 , z)/ip{z)) Im(v9(^2, z)/ip{z)) 



lvsnp{z) 



Im(l/^(z)) 



By Lemma 4.1, if{(^i,z)/ip{z) and Lp{(^2-, z) / ip{z) are Stieltjes functions, and by Proposi- 
2.1[ c), \/lp{z) is a complete Bernstein function. Hence, the second summand on the 



tion 



right-hand side of (4.1) is non-negative when Im^ > (Proposition 2.1[ a,b)). The first 
one is non-positive, but below we prove that it is bounded below by Im(y9(^i, ^2, z). 



For z = X + iy with real x, ?/, we have as in the proof of Lemma 4.1 



Im V3(^i, z) Im (y9(^2, z) - Im ip{z) Im V3(^i, 6, z) 



00 fOO 



^0 



X 



e^ + ci vei + ci a + ci 



((vP(Cl)+x)2 + y2)((vI;((2)+x)2+y2; 



dC,idC,2 
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By a similar symmetrisation procedure as in the proof of Lemma 4J^, we obtain that 
2 Im (^(^1, z) Im ip{^2, z)-2lm (p{z) Im (p{^i, ^2, z) 

-2 

52 



1 fca noo / f2 f2 \ / f2 f2 

-"-/// SI SI \ / S2 S2 



X 



((M/(Cl)+x)2 + y2)((vI;(^2)+3;)2 + ^2) 

and again the integrand on the right-hand side is non-positive. Hence, 

lmip{^i,z)lmip{^2,z) -lmip{z)lm(p{^i,$,2,z) < 0. 

When y > 0, we have lmip{z) < 0, and so 

lm(p{^i,z)lmip{^2,z) 

— > lmip{^i,^2,z). 

Lm(p[z) 

By ( 4.1[ ), it follows that lm(^(^i,^2, z) < whenever Imz > 0. Clearly, Im(^(^i, ^2, 2;) = 
— lm(^(^i,^2,z) > when Im^; < 0. By Proposition |2.lK b), it remains to prove that 
^2, z) >0 when z > 0. 
Since ip(^i,C,2,z) is real when z > and lmip{^i,^2, z) < when Imz > 0, we have 
id / dz)(p{^i,^2i z) < for z > 0. It follows that ip{^i,E,2iZ) is a non-increasing function of 
2; > 0, and hence it suffices to show that lim2_^oo '^(G? ^2, z) > 0. When z — )■ 00, by the 
definition of ip{^,z) and 99(^1,^252;) and dominated convergence, the functions ip[E,i-,z) 
and ip{^i,^2^z) converge to 0. Also, ip{^2^ z) / ip{z) < 1 for all z > 0. We conclude that 
(^(^1, ^2, z) converges to as z — )■ 00, and the proof is complete. □ 

For A, ^ > 0, we define 

1 f2 ^/MN p2 



and, as in (3.3), 



1 p ^'{\)di 

Kx = pv ' 



7r 







We denote the boundary limits of f{z), (p{C,, z) and ip{C,i, ^2, z) along (—00, 0) approached 
from the upper half-plane by </)+(— z), ip~^{^, —z) and ¥^^(^1,^2, —z) {z > 0), respectively. 
The existence of these limits is a part of the next result. 

Lemma 4.3. For ^, A > 0, we have 

^ ^'^^^^^ = ^ ^ — n>^) — ■ 

Furthermore ( cf. Lemma 4-1^ , 

(^+(e,-^(A)) KxLxiO-KxiO 



Im 



and for ^1,^2, ^ > (cf. Lemma 4-. 2), 



Im<^+(ei,6,-^(A)) 

1 KxLxiii) - Kxi^i) KxLx{^2) - Kx{^2) 



^'(A) ^/iTkI vTTlr 
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Proof. For A > and any bounded, continuously differentiable function / on (0, oo), we 
have (see Section 1.8 in |18| ) 

/•oo ^ 
lim / T77^ 77T^ ^ fiC)dC 



/"°° 1 z 

Hence, the boundary hmits \&(A)), f^iC,, — \I/(A)) and ^2, ~^('^)) (approached 

from the upper half-plane; here and below ^, ^i, ^2 > 0) are continuous functions of A > 0, 
given by 

y,+(-v^(A))=--pvr— 1— ^ ' 



1 /-oo f2 1 • f2 

¥'ne,-^(A)) = — pv/ 7^ ;^77TX-^^;77T ^ 



vr'io e + e*(A)-v&(C) ' v[/'(A)e + A2 
_ irA(0-zLA(0 
^'(A) 

and (we only need the imaginary part) 

Im^ (ei,6,-^(A)) = -^^2^^2^ = 

We obtain that 

(^+(e,-^(A)) i^'AL,(0 - i^A(e) 



Im 



(^+(-vl/(A)) 1 + Kl 



and therefore, by (4.1), 



V.+ (ei,-vl^(A))y.+ (6,-^(A)) 
V.+(-vl/(A)) 

^A(ei)^A(6) , i^A^A(ei) - Kxi^l) KxLx{i2) - i^A(6) 1 + Kl 



^'(A) l + ^I 1 + ^! *'(A) 

1 irAi^A(ei) - Kxii,) KxLx{i2) - i^A(6) _ ^^(^^)^^(^^) 



^'(A) \^ ^A+^ ^/^+Kl 
The proof is complete. □ 
Lemma 4.4. For a// ^1,^2 > and t > 0, let 

^ r A-.L.(6)-A-.(6) A-.L.(fe)-A-.(fa) „, 

Then 

^(6,6,0) = ^-^, (4.3) 



anc? /or all z > 0, 



1 Z""" 

- / $(6,6,i)e^^*cit = <^(ei,6,^)- 

71" Jo 



(4.4) 
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Proof. Note that by Lemma 4.3, we have 



Im<^+(ei,6,-^(A))e-**(^)^'(A)t/A 
Im(^+(^i,^2,-s)e"*'(is. 







By Lemma 4.2, lm(p^(^i, —s) ^ 0. Hence, by Fubini 



POO 

Jo 



s + z 



Im<^+(^i,6, -s)ds. 



By the definition, the coefficients Ci,C2 in the representation (2.2) of the Stieltjes func- 
tion (p{^i,^2,z) vanish. Hence, by Proposition 2.3K b) and the inequahty (^(^1,^25-2^) < 
f{^iy^2,z) {z > 0), the coefficients Ci,C2 in the representation (2.2) for ip{C,i,C,2,z) vanish 
too. By (2.2) and another apphcation of Proposition |2.3[b). 



s + z 



Im^-^ {^i, ^2, -s)ds = Tnp{^i,^2,z) 



This proves (4.4). 

By monotone convergence, $(^1,^2?^) is right-continuous at t = 0. This and (4.4) 
imply that 

/■oo 

$(^1,6,0) =\imz ^^i,^2,t)e-''dt = 7T \im{z^{^,,^2, z)). 

As z —> 00, by monotone convergence, the function z(p{^, z) converges to ^, but z(p{z) 
diverges to infinity. Hence, z(f{C,i,z)(f{C,2,z)/(f{z) converges to 0. It follows that (see the 
definition of (^(^1,^2,-2) in Lemma 4.2[ ) 



$(6,6,0) = TT lim(z(^(6,6,2;)) 



monotone convergence was used again for the second equality. The integral can be easily 
evaluated (we omit the details), and (4.3) follows. □ 



With the above technical background, we can prove generalized eigenf unc tion expansion 
of the operators P/^^*^. Let Heven, ^odd and H be defined as in Theorem 
prove that the operator Heven is isometric, up to a constant factor. 



L9 



First, we 



Lemma 4.5. For any even functions f,gE L^(R \ 0), 

(Heven/, nevcnfi') ^2 ((0,00)) = ^ (/, d) L^{n\0) " 

Proof. Let eg(x) = e~^'^', x G R \ 0, ^ > 0. By Theorem L4, for 6 A > we have 



Hevene5(A) 



FA(x)e-«l^lda; 
-2KxiO cos^a + 2La(0 sin^, 

2 K,L,iO - K,iO 



(4.5) 
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By (|42|) and (Q, it follows that for 6,6 > 0, 

4$(6,6,0) 



even 



66 

2n 



L2(R\0) • 



6 + 6 

The family of functions eg is linearly dense in the space of even L^(R \ 0) functions. The 
result follows by approximation. □ 

The following side-result is interesting. It is obvious when "^l^) = for a complete 

Bernstein function ip. In the general case, however, direct approach seems problematic. 



Proposition 4.6. We have 



Fx{x)e-^\''\dx > 0, 



A,e>0. 



Proof. By Lemma 



4.3 



and (4.5), we see that for e^{x) = e' 



Vr+g v^+(6-^(A)) 



which is nonnegative by Lemma 4.1 



□ 



The space L^(R\0) is the direct sum of the spaces Lg^gj^(R\0) and L^^j(R\0), 
consisting of even and odd L^(R \ 0) functions, respectively. Recall that Ilodd is the 'sine 
part' of the Fourier transform. Hence, 7r~^/^nodd is a unitary mapping of L^j^(R \ 0) onto 
L^((0, c>o)), and Ilodd is a zero operator on ^^^^^^(R \ 0). In a similar manner, 7r~^/^neven is 
an isometric mapping of L^^^^ill \ 0) into L'^{{0, oo)), and Ilgven vanishes on L^^j(R \ 0). 
It remains to prove that 7r~^/^neven is a unitary operator on ^^^^^^(R \ 0) onto L'^{{0, oo)). 
An equivalent condition is that the kernel of the adjoint operator 11*^^^^ : L^((0, oo)) — )■ 
L^(R \ 0) is trivial. This is proved in the following result. 



Lemma 4.7. For all f,g e L^{{0, oo)), 

(n 

even/ 5 -^evens') (R\0) 

and 



L2((0,oo)) ' 



t ■'■-'■even 



/ 



n* (e 

even \^ 



/), 



t > 0. 



(4.6) 



(4.7) 



Proof. As in the proof of Lemma 5.2 in [12j, formula (4.7) for / G Cc((0,oo)) follows 
directly from Theorem L4 and Fubini. Indeed, let supp/ C [Ai,A2], < Ai < A2. 
Using (1.2), (1.4) and Fourier inversion formula, it is easy to prove that F\{y) is bounded 

uniformly in A G [Ai, A2] and ?/ G R \ (we omit the details). Furthermore, (x, y) is 
integrable in ?/ G R \ 0. Hence, 

{x,y)Fx{y)f{\)d\dx 



t ■'■■'■even 



/ 



R\0, 
Pt ' 
R\0 ^Ai 

'%-**WF,(x)/(A)rfA 

Ai 



H! 



Since n*^g^ is a bounded operator, (4.7) extends to general / G -^^^((0, 00)) by an approx- 
imation argument. 
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The adjoint of an isometric operator is isometric if and only if the kernel of the adjoint 
is trivial. Hence, by Lemma 4.5, it suffices to show that the kernel of Hl^en trivial. 

Suppose that / G L'^{{0, oo)) and 11*^^^^/ = 0. We claim that this implies that / = 0. 
By n*^^^(e-**/) = for any t > 0. Let e^{x) = e"?!^!, x eR\0. We obtain that 



for alH > and ^ > 0, 

= (n* 

" \ evf 



L2(R\0) 



L2((0,oo)) 



By substituting s = ^I^(A), we obtain that 



f{^-\s))U,.,^e^{^-\s)) 
e as 



-t*(A) 



0, 



/(A)n 

even 



,0 vl/'(vl/-i(s)) 

It follows that for each ^ > 0, /(^"^(s))nevene^(^"^(s)) = for almost all s > 0. 

In other words, for each ^ > 0, f {X)IlevenG^{X) = for almost all A > 0. By Fubini, 
for almost every A > 0, /(A)nevene5(A) = for almost all ^ > 0. But for each fixed 
A > 0, neveneg(A) (as a function of ^ > 0) is the Laplace transform of a nonzero function 
(namely, 2F\(x)l(o,oo)(a;)), and so it cannot vanish for almost all > 0. This implies that 
/(A) = for almost every A > 0, proving our claim. □ 



Proof of Theorem 1.9. The operator 11/ = (Ileven/; nodd/) is unitary by Lemmas 4.5 



and 4.7, and the discussion preceding the statement of Lemma [4. 7[ 

Let / G L^(R \ 0), and let / = /even + /odd be the decomposition into the even and odd 
part. Clearly, P^'^^feven is an even function, and P^'^^fodd is an odd function^ Hence, 

Heven/odd = H 

even Pt /odd = and Hodd/e 



Hodd-Pt'^^^/e- 



0. By Lemma 



4.7 



H, 



I-' t jeven 



TT-^He 



n* n f 

even even J even ; 



TT 



(U* fe"**H f )) = e""^H f 

V evenV^ -'--'-even J even/ / ±±QYenJ • 



(4.8) 



Furthermore, by the strong Markov property, pt{x — y) — Pt {x, y) = Ex{pt{y)] tq < t) 
(the Hunt's formula), and the right-hand side is an even function of y (for any x G R \ 0) 
Therefore P^'^'^fodd = Ptfodd- By the remark preceding the statement of Theorem 



see also Remark 2.7), it follows that 



1.9 



nodd-Pi^^°/odd — Hodd-Pt/odd — 6 *'''Hodd/odd 



The above properties combined together prove (Lll) 



Proof of Corollary [T^ Let / G C^(R\0). Then H/ G L'^{{0,oo)) 
hence e~*'^H/ is a pair of integrable functions. Hence, 

prm 



(4.9) 

□ 

L2((0,oo)), and 



7r-iH(e-**H7)(a; 

%-**Wf,(x) f{y)F,{y)dy]dX 
1 



+ 



TT 



-t*(A) 



sin(Aa;) 



fiy) sin(A?/)(i?/ dX 



By Theorem L8[ a), Fx{y) is bounded uniformly in A > and y G R \ 0. By Fubini, it 
follows that 

Pl^'^'fix) = - / fiy) / e-**(^)(FA(x)FA(i/) + sin(Ax) smiXy))dXdy. 

J-oo Jo 
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TT 



e-**(^) sin(Ax) sm{\y)dX 



"**(^)(cos(A(x - y)) - cos(A(x + y)))dX 



-oo 
oo 



1 
1 

Ptix - y) -ptix + y). 



-t^(X)( -i\{x-y) __ -iX{x+y) 



)dX 



□ 



Proof of Theorem 1.11. Recall that we need to prove that 



^ 1 P cos^9Ae-**(^)^'(A)FA(x) „ 

TT Jq m{x) 



(4.10) 



Note that Gx(fS) = sin-i^A — -^a(O) = sin-t^A- For A, t > and a; G R \ we decompose the 



integrand in (4.10) 



cos i?Ae-**(^)^'(A)FA(a;) e-**W^'(A) 



M/(A) 



M/(A) 



^sin(|Aa;| + -(^a) — sin-t^A) 



+ ^P^(Ga(0)-Ga(x)). 



(4.11) 



By the assumption, \l/(0 = where ip is non-negative, increasing and concave. 

Hence, ^'(A) = 2Xij\X'^) < 2X{ijlx'^)/ X'^) = 2^(A)/A. It follows that 



e-**W^'(A) 



M/(A) 



^sin(|Aa;| + i9a) — sin-d). 



< 2|x|e-*'"(^). 



Furthermore, the second summand on the right-hand side of (4.11) is non-negative by 



Theorem 1.8 It follows that the integral in (4.10) is well-defined (possibly infinite). 



Denote the right-hand side of ( |4.10[ ) by f{t,x) (t > 0, x G R \ 0). If eg(x) = e'^l""', then 
for ^, 2; > we have, by Fubini, 



/(t,x)e-«l^le-"*dxdt 



1 /'°°COS^A^'(A)ne.ene^(A) 



J -00 



{z + ^{X))^{X) 



dX: 



the use of Fubini here is justified by the decomposition of the integrand into an absolutely 



integrable part and a non-negative part, as in (4.11). Using cos'i^A = (1 + K^) and 
formula (4.5), we obtain that 



/•oo poo 

/ / f{t,x)e-^\''\e-'*dxdt 

Jo J -00 
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By Lemma 4.3 substitution \I/(A) = s and Proposition 2.3[ b), 



r( ^ Im^^:^ii^^vI;YAWA 



2 ^li^ Vii.e) vii^z) 



^z \s^o+ (p{e) (p{z) 

Note that ip{C,,z) < (p{z), so that the above expression is finite. In particular, f{t,x) is 
finite for almost alH>0,a;GR\0. 

On the other hand, for ^ > and x G R, denote gz{x) = E^,. exp(— 2;ro). Then 
gz{x) = Uz{x)/uz{0), where Uz{x) is the resolvent (or 2;-potential) kernel for the operator 
A (see Preliminaries). Hence, 



where, by the Fourier inversion formula 



Furthermore, for z > and x G R \ 0, 

E^exp{-zTo) gz{x) 



Jo 

By Fubini and Plancherel's theorem. 



POO POO 

/ / Px(ro < t)e'^^'=\e~''dxdt 

Jo J -oo 

-I l*00 r\ POO ^ 

= - / gzix)e-^^^\dx = — / ^^Tg^iQdC 

Z J-oo Jo i + 



Txiz Lp{z) Jo e + C ^(C) +Z ^ iz ^{Z 

In particular, we have 

/•oo POO 

Jo J-oo 

Z J —oo 

= lim y ^ P-(^o < t)e-^^^^e-''dxd?j 



2 v^(e,^: 



£^0+ v9(£:) 
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and so finally 



/•oo poo 

/ / P^(t < To < oo)e-«l^le-"Vxrft 

Jo J -oo 



lim 



f{t,x)e-^\''\e-'*dxdt. 



By the uniqueness of the Laplace transform, formula (4.10) follows for almost every pair 
t > 0, x G R \ 0. By dominated convergence, both sides of (4.10) are continuous in t > 0, 
and the theorem is proved. □ 



Remark 4.8. One could try the following argument for the proof of Theorem |1.11[ By 
monotone convergence, 

{x,y)dy 



Px(ro>t)= I 

J —c 



R\0/ 
Pt ' 



lim 



/oo 
^^■^^(a;, |/)e~^'^'(i?/ = lim Pf^^°e^(a;) 
-oo «^0+ 



where e^{x) = e ^'"^'L For each ^ > 0, by Theorem 



1.9 



5-s>0+ 

we have 



,R\0 



e^[x) 



n* (e"**n 



,eg)(x) 



liX)d\. 



(4.12) 



Recall that as ^ — )■ O"*", Ueven^^i^) converges to cos'i9A^'(A)/\I'(A). Hence, if the above 
integrals were uniformly integrable as ^ — )■ 0+, we would obtain 



Px(ro > t) 



cos^9a*'(A) 



d\. 



However, if Pxi^o < oo) < 1 (that is, when Xt is transient), we know by Theorem 1.11 



that this formula is false! This proves that the integrand in (4.12) in some cases fails to 
be uniformly integrable, and its limit (in the sense of distributions) may contain a point 
mass at 0. For this reason, a more direct approach to the proof of Theorem 1.11 , similar 
to the argument sketched above, seems to be problematic and require stronger regularity 
conditions on □ 



5. Examples 



In this section we consider a few examples: processes frequently found in literature and 



an irregular case. We focus on spectral theory, that is, on Theorems 1.4, 1.8 and 1.9, and 
only record formulae given in Corollary 1.10 and 1.11 As mentioned in the introduction. 



detailed properties of first hitting times to points will be studied in a separate article. 

Plots have been prepared using Mathematica 8. 1 at the Wroclaw University of Technol- 
ogy. Numerical integration was used for the computation of d\ and fast Fourier transform 
for the calculation of Gx{x). 

Example 5.1. Let ^E'(^) = |^|", a G (1,2], be the Levy-Khintchine exponent of the one- 
dimensional symmetric a-stable Levy process. The distribution of the hitting times to 
points for these processes have been previously studied in 
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'&\{a) = ir/a — tt/2 



Fi{x) (blue), Gi{x) (purple) and 




sin ( a: 


+ i?i) (brown) for a = 1.1 


/ \o.75 




/ V'' 




/ 








/ -0.25 




/ ""-^ 




-0.75 
-1 





Fi{x) (blue), Gi{x) (purple) and 
sin(|x| + iDi) (brown) for a — 1.5 



Fi{x) (blue), Gi{x) (purple) and 
sin(|a;| + t?i) (brown) for a — 1.9 



Figure 1. Plots for 'i'(^) ^ |^|". 
By substituting ^ = As in the integral over (0, A) and ^ = X/s in the integral over 

1 



(A,oo), we obtain (see (3.3)) 



Kx = -- pv 



^ Jo A" - e 

By series expansion and Fubini, 

OO / /"I 



a 

TT 



1 



1 - - S") 



ds. 



n=0 ^ ' ra=0 ^ ^ ' 

Recall that for z G (0, vr), 

7rcot(7rz) = — h > 1 . 

z ^-^ \ z — n z + n I 

n=l ^ ' 



Hence, after some rearrangement of the sum, 



Kx 



vr ^ V 1/a ■ 



n=0 



1/a + n — (n + 1) 



+ 



Hence, 



= arctan(KA) 



— cot 



vr 



vr vr 
a ^ 2 ' 



(5.1) 
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and cos'i^A = sm(7r/a), sin(^9A) = — cos(7r/Q;). Furthermore, by (1.9), 



^ , , aA" ^ sin(7ra/2) sin(7r/a) 
Gxix) = ^ — '-^ ^^^x 



TT 



X 



a sin(7ra) sin(7r/Q;) 



X2a _ 2X^^» cos(7r«/2) + 



Hence, 



F\{x) = sin ( |Aa;| + 



TT 



vr TT 



1 — 23°" cosfvra) + s 



2a 



a 



a sin(7ra/2) sin(7r/Q!) 



X 



TT 



X 



1 - 2s" cos(7ra/2) + s 



2a 



(5.2) 



Note that Fx{x) = Fi{Xx), which could be deduced a priori using scahng properties of Xj. 
As expected, in the hmiting case a — )■ 1+ we have -t^x — )■ tt/2 and Gx{x) — )■ (the Cauchy 
process never hits 0, hence Fx{x) = cos(Ax) = sin(|Ax| + tt/2) is the even eigenfunction 
in R \ 0). Also, for a = 2, we obtain = and Gx{x) = {Fx{x) = sin |Ax| is the even 
eigenfunction for the Brownian motion in R\ 0). See Figure [T] for plots of Fx- 

Since Xt is recurrent, we have Px{^o = oo) = 0. Hence, by Theorem 1.11 and scaling. 



Px(ro > t) 



a sin(7r/a) 



TT 



oo Q-tX" 



FiiXx) 



d\. 



By Corollary 1.10 



e~*^"(Fi(Ax)Fi(A?/) + sin(Aa;) sm{Xy))dX. 



These formulae are applicable for numerical computation. 

Example 5.2. Let a G (0,2), /3 > 0, and ^(^ = + The corresponding Levy 

process Xf is a mixture of Brownian motion and the symmetric a-stable process. Note 
that Xt is recurrent if and only if a > 1. We have 



'dx = arctan ( pv 



2A + a/3 A 



a-l 



A2 + /3A" - e - /3e 

This integral does not have a closed form, but it can be proved that limA-s.oo ''^a = and 
limA^o+ i^x 

Fxix) = 



mm{Tc/a — tt/2,tt/2). By Theorem 1.8 
(2A + a/3A"~i)cost9A 



sin(|Aa;| + -^x 

POO 

X 



TT 



sin(a7r/2)^^ 



(A2 + /3A" - ey + - 2 cos(a7r/2)(A2 + /3A" - e)^' 



dC 



Example 5.3. Suppose that \E'(^) = (^^ + 1)"^^ — 1 is the Levy-Khintchine exponent of 
the relativistic a-stable Levy process, a G (1,2), corresponding to unit mass. In this case 
again 



arctan 



pv 

TT Jo 



a 



A(A2 + i)"/2-i 



(A2 + l)°/2 - (^2 + 



d^ 
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1 2 3 4 5 6 



(blue) and (purple) 




Fi{x) (blue), Gi{x) (purple) and 
sindxl (brown) 



^4(0;) (blue), 6*4(2;) (purple) and 
sin(4|a;| + 1?4) (brown) 




-^4.5(2;) (blue), 6*4.5 (x) (purple) and 
sin(4.5|a;| + 194.5) (brown) 



F^^x) (blue), G^^x) (purple) and 
sin(5|x| +1^5) (brown) 



Figure 2. Plots for ^'($) = ,^2_^9(l-cos^). Note the differences between A = 1 (i?a > 0, 
Gx attains global maximum at 0), A = 4 (i9a > 0, G\ has a positive local minimum 
at 0, A = 4.5 (i9a and G\{0) close to 0, slow convergence of Fx{x) to sin(|A2:| + dx) as 
\x\ 00) and A = 5 (i^a < 0, Gx has a negative global minimum at 0). 



cannot be evaluated symbolically. However, with some effort, one can show that 
limA^.o+ "i^A = and limA^-oo i!^x = n/a — 7^/2. We also have 

^ , , , , , «A(A2 + l)"/2-icosi?A 
Fx{x) = sin(|Ax| + ^^a) ^ 

X r sin(W2)(l-a°/^ 

A (A2 + l)" + (l-e2)"-2cos(«7r/2)(A2 + l)-/2(l-^2)a/2'^ 

In particular, -Fa(x) — sin(|Ax| + "^x) decays exponentially fast as |x| — )■ 00. 
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Example 5.4. Let = + 9(1 — cos^). The corresponding process Xt is the sum of 
the Brownian motion (with variance 2t) and an independent compound Poisson process 
with jumps ±1 occurring at rate 9. It can be easily checked that \£''(^) > for ^ > 0, 
and therefore Theorem 1.4 applies. However, -^x is negative for some A > 0, and G\{x) 
is a signed function, see Figure [2j 



Since = \t'(v^) = ^ + 9(1 — cos v^) is not concave. Theorem 1.11 (or Theorem 1.1 ) 
cannot be applied to Xt. 



Example 5.5. Consider the truncated symmetric a-stable Levy process Xt, that is, a 

for some c > 0. Then, 

1 - C0S{^/^X) 12 f"^ 1 - COS S 



pure-jump Levy process Xt with Levy measure c\x\ ^ "l(_i^i)(x) for some c > 0. Then, 



J-1 l-^l Jo 



ds. 



Clearly, > for all ^ > 0. We claim that ^/'"(O < for all ^ > 0. We have 

2V'"(0 sinVe ^ l-cosy/^ a{2-a) f^^ 1-coss 
^^ = ^37^-(2-a) io s^+- 

Since sins < s — + s^/120 and 1 — coss > — s^/24 for s > 0, we have 
2<(0 ^1 2-a^2-a a ^ a{2 - a) 



^ 6 120 2^ 24 2^ 24(4 -a) 
1 e 



3(4 - a) 120 

and hence V'"(0 — ^ — 40/(4 — a). When ^ > 40/(4 — a), then in particular 
^ > \/l2, and therefore, by ( |5.3 ), 

2V^^X0 «(2-a) /"^ sV2 - s724 



^ 1 12i-"/2c 



^3/2 ^ (4 - a) ^ y - ^3/2 

Our claim is proved. 

In particular, we may apply Theorems 1.4, 1.9 and 1.11 to truncated stable processes. 
Due to a rather involved description of \E', the computations are significantly harder in 
this case. See Figure |3] for plots of 'dx and Fx. 

Acknowledgements. I would like to thank Kouji Yano for suggesting me the study of 
the spectral theory on R \ and interesting discussions on this subject. I also thank 
Jacek Malecki for inventing and explaining me the proof of Theorem 1.9 in ^3j, which 
made the present proof of Theorem L9 possible. 
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dx {a =1.1) 



z?A (a = 1.5) 




^4(0;) (blue), Gi{x) (purple) and 
sin(4|a;| + 1)4) (brown) for a = 1.1 




Gx{x) for a = 1.1 and A = 1/16 
(blue), A = 1/4 (purple), A = 1 
(brown) and A = 4 (green) 
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